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We examine one of the advantages of Ashtekar's formula- 
tion of general relativity: a tractability of degenerate points 
from the point of view of following the dynamics of classical 
spacetime. Assuming that all dynamical variables are finite, 
we conclude that an essential trick for such a continuous evo- 
lution is in complexifying variables. In order to restrict the 
complex region locally, we propose some 'reality recovering' 
conditions on spacetime. Using a degenerate solution derived 
by pull-back technique, and integrating the dynamical equa- 
tions numerically, we show that this idea works in an actual 
dynamical problem. We also discuss some features of these 
applications. 

PACS No.: 04.20. Cv, 02.40.-k, 04.20.Fy, 04.25.Dm 



I. INTRODUCTION 

A decade has passed since the proposal of the new for- 
mulation of general relativity by Ashtekar [El. By using 
this special pair of variables, the framework has many 
advantages in the treatment of gravity. The constraint 
equations which appear in the theory become low-order 
polynomials, and the theory has the correct form for 
gauge theoretical features. These suggest possibilities for 
treating a quantum description of gravity nonperturba- 
tively. 

Here we examine another advantage of the SO (3)- 
ADM (Ashtekar) formulation of general relativity: the 
ability to dynamically evolve a spacetime with a degen- 
erate metric consistently. A 'degenerate point', here, is 
defined as a point in the spacetime where the density e 
of 3-space vanishes, as we denote in §[n| This advan- 
tage comes from the fact that all equations in Ashtekar's 
formulation do not contain the inverses of the variables. 

There are several motivations for studying 'degenerate 
points'. The first one comes from a quantum cosmo- 
logical description of the very early history of the Uni- 



verse. One scenario describes that our Universe was born 
and evolved in Euclidean spacetime, then emerged with 
Lorentzian metric through a quantum tunneling process. 
Kodama Q found an exact solution for Bianchi type IX 
spacetime using Ashtekar's framework, which has both 
Euclidean and Lorentzian sections connected by an ana- 
lytic continuation. Apparently, such a signature changing 
process should pass through a degenerate point, but so 
far the dynamics of this process remains unclear. 

A second motivation comes from the gauge theoretical 
point of view. We can set a finite gauge transformation 
so as to have a degenerate point of spacetime from a non- 
degenerate one naturally |3). Therefore when we calcu- 
late the functional integrals, we have to include degener- 
ate solutions due to the gauge invariance, and studies on 
degenerate points will be indispensable to construct the 
complete theory of quantum gravity in the near future. 

A third interesting application exists in topological 
field theory. Since Ashtekar's dynamical variables E l a 
and Af are assigned positive mass dimension one, a short- 
distance limit (i.e., region of quantum gravity) must have 
degenerate phase: E = A = 0. Such a degenerate phase 
is called the topological (unbroken) phase of quantum 
gravity, in which diffeomorphism is unbroken M. Find- 
ing how to evolve from the degenerate unbroken phase 
to the non-degenerate phase is one of the important is- 
sues, and it is expected to be a key to construct quantum 
gravity in four dimensions. 

Another possibility which makes degenerate points in- 
teresting is the evolution of classical (Lorentzian) space- 
time. As already pointed out by some authors (|||, the 
Ashtekar formulation is also attractive for numerical rel- 
ativity. If this tractability of a degenerate point works 
also in this context, then we will be able to analyze focus- 
ing or shell-crossing features or (coordinate) singularities 
in spacetime. 

In this paper, we study what problems appear and how 
to solve them when we treat the dynamics through a de- 
generate point. Here, the dynamics means the evolution 
of the three- hypersurfaces expressed by Ashtekar's vari- 
ables (E % a , Af), fixing the gauge freedom of the lapse func- 
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tion N, shift vector iV* and triad lapse function Aq (this 

name appears in j5j) every time step. This is a common 
procedure for treating the dynamics in the ADM formu- 
lation, though we have an additional gauge freedom. In 
order to set up a problem of dynamics clearly, we fix our 
spacetime to have Lorentzian signature. 

In the last few years, several studies have been done 
on degenerate metrics. Some solutions have been found 
by connecting degenerate and non-degenerate metrics J7j- 
P], and degenerate examples in terms of Ashtekar vari- 
ables are shown in [ |To[ . Some alternative approaches over 
the framework of Ashtekar were considered in and 
JL2| . However there are few comments on the dynamical 
evolution of a degenerate metric. Bombelli and Torrence 
|l3f commented on the conditions for 'passing' degenerate 
points. They speculate the degeneracy of lapse function, 
condition for finite 'passing' coordinate-time, and possi- 
bility of divergence of variable s at t he degenerate point. 
We will develop their idea in §111 A. 



Assuming that both all dynamical variables and the co- 
ordinate time are finite (passing condition in §[n]) , we seek 
the condition which enables us to pass degenerate points 
for the case of a vacuum spacetime with/ without a cos- 
mological constant. We take mainly two approaches. The 
first approach, we named intersecting slice approach, pur- 
poses that a foliation passes a degenerate point directly. 
This is a natural passing behavior which we expect from 
a context of Ashtekar's formulation. However, we show 
that such a direct treatment does not work normally both 
in ADM and Ashtekar formulation, whether we impose 
the reality conditions or not. This is described in §111. 
The second approach, we call deformed slice approach, 
takes a complex path around a real degenerate point. 
This approach is described in § |Vl in which we propose 
some conditions for recovering the real manifold after a 
detour around a degenerate point by complcxification. 
In §[v|, we derive a degenerate solution by pull-back tech- 
nique and show that the deformed slice approach works 
well for 'passing' a degenerate point using a numerical 
integration of the dynamical equations. We conclude 
that the essential trick for 'passing' a degenerate point 
in Ashtekar's formulation exists in their complexification 
of the basic dynamical variables. We devote §V1 to dis- 
cussion. 

When we apply the Ashtekar formalism in classical 
general relativity, we need to impose reality conditions. 
We treat these reality conditions based on work devel- 
oped by two of us Q, and some essential points are at- 
tached in the Appendix along with a full description of 
the notations we use throughout this paper. 

We use greek letters (fx, v, p, ■ ■ •), which range over 
the four spacetime coordinates 0, • • • , 3, while uppercase 
latin letters from the middle of the alphabet (I, J, K, ■ ■ •) 
range over the four internal SO(l,3) indices (0), • • • , (3). 
Lower case latin indices from the middle of the alphabet 
(i,j,k,...) range over the three spatial indices l,---,3, 
while lower case latin indices from the beginning of 



the alphabet (a, b, c, ...) range over the three internal 
SO(3) indices (1), • • • , (3)0. We use volume forms e a bc'i 



^abc c>| 



II. DEFINITIONS OF 'PASSING' 

In the beginning, let us clarify some terminology. The 
basic variables and equations in Ashtekar's formulation 
are summarized in the Appendix. 

A 'degenerate point', in this paper, is the point in the 
spacetime where the density e of 3-space vanishes (not 
4-space). In the Ashtekar formulation, the density is de- 



fined as e = y det E\ , which corresponds to the same 

condition in the ADM formulation as e 2 = det 7^ , where 
7ij is 3-metric. We ask readers to remember that such 
a degenerate point does not always mean a physical sin- 
gula ri ty. In the basic equations neither t he c o nstr aints 
(|A5|)-( A7) nor the dynamical equations ( A8)-([A9|) in- 
clude any inverse of the dynamical variables (Af and E % a ), 
even if we include a cosmological constant g. This fact 
suggests to us that we can 'pass' such a degenerate point. 
That is, we expect that all calculations can be continued 
even if we have a degenerate metric during the time evo- 
lution. 

Note that, in the ADM formulation, this is impossible 
because the equations in the ADM include an inverse of 
the variables. We also remark that we can not transform 
the Ashtekar variables onto the ADM three-hypersurface 
£ at the degenerate point. This is because the 3-metric 
7y is given by 7^ :— e°e°, where e.f :— (-E^) -1 and 
Eq := E l a /e, and the last quantity diverges as e — > 0. 
Therefore, at the degenerate points, one can not recover 
the ADM variables from Ashtekar variables. 

In order to say 'pass' degenerate points, we require the 
following four passing conditions: 

(a) Ashtekar variables E l a ,A°;, N, N 1 ,Aq must remain fi- 
nite throughout the calculation, 

(b) the spatial derivatives of them must also be finite 
throughout the calculation (because they appear in the 
equations of motion) , 



(c) all the constraints ( A5),( |A6|) ,( A7) and the equations 
of motion ©,@ are satisfied, and 



1 We raise and lower fi, v, p by g^ v and (Lorentzian met- 
ric); I,J,K by n IJ = diag(-l, 1, 1, 1) and rju; i,j,k by j lJ 
and 7ij (3-metric). 

2 The cosmological constant term in (A8) looks as if it might 



be divergent, but the relation eef = (l/2)e abc € ijk E J b E^ guar- 
antees its finiteness. 
3 From the assumptions (b) and (c), we have that the time 

differentiation of dynamical variables E % a , must also be 
finite. 
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(d) the calculation must be finished in finite coordinate 
time. 

The problem is whether we can make dynamical fo- 
liation through a degenerate point under these passing 
conditions (a)-(d) above. In the following sections, we 
take two approaches. The first one, which we call the 
'intersecting slice approach', attempts to pass a degen- 
erate point directly, and the second one, 'deformed slice 
approach', takes a foliation in complex region. 



III. INTERSECTING SLICE APPROACHES 

In this section, we describe the possibility of an 'inter- 
secting slice approach', which attempts to pass a degener- 
ate point directly. Here 'direct' means that the dynamical 
Ashtekar variables run into the degenerate points, never 
detour the points. We assume that N and N 1 are real as 
usual. To say our conclusion first, we face at least two 
problems, which we call "cusped lapse/density problem" 
and "divergence problem" , the latter requires severe con- 
ditions. 



A. Cusped lapse/density problem 



call this the cusped lapse problem. We note that this is 
not a serious problem to execute, because the choice of 
lapse in such a cusped form is only within a freedom of 
gauge, although there remains its natu raln ess of foliation. 

In terms of t, however, the request (3.1) is not restric- 
tive as follows. When N = dr/dt = r s (0 < s < 1), we 
see t = r 1 - s /(l - s), thus N = (s + 1)^/(1-^/(1-*). 
The power s/(l — s) can be taken arbitrary positive even 
if we require < s < 1. The only request here is N = 
at t = 0. So the lapse is not necessarily cusped in terms 
oft. 

Furthermore we will see that the density |e| should 
also satisfy a condition which is similar to the cusped 
lapse condition. We call this problem the cusped density 
problem. We assume e = 0att = r = and \N\ < M 

(bounded). Since N > 0, we have eN = \eN\ = |e| \N\. 
Then we see 



00 > 



eh 
~N 



dr 



eh 



> 
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-r eN J- To \e\ \N\ M J_ To \e\ 



dr 



Similarly to the lapse, the density is restricted in terms 
of t as |e| = t s (0 < s < 1), and is not restricted in terms 
oft. 



A troublesome variable in an evolution through a de- 
generate point is the (inverse) densitized lapse N := 

N/es, where = «/ det jij ■ Since N is held finite [con- 
dition (a)], the ADM lapse N = Ne^ vanishes at the 
degenerate point es = 0. Notice that the ratio of the 
proper time r to the coordinate time t is not N but N. 

Thus we are afraid that the calculation exhausts an in- 
finite amount of time, i.e., condition (d) is violated. Let 
us take t and r such that r = t — at the degenerate 
point es = 0. Then condition (d) is denoted by n 



St 



dT 



dt= - < 00. 

J -TO N(T) 



(3.1) 



Here the range of the integral is arbitrary but includes 
the zero point t = r = 0. Note that iV > since it is 
a lapse function. This condition originally appeared in 

In terms of t, the request ( |3.1| ) is considerably restric- 
tive because N — at r = 0. For example, let us consider 
the form N — \t\ s where s is constant. The condition for 
(O) is < s < 1, e.g., N 



\t\ satisfies fl3.l|) . We see 
dN/dr does not exist at r = This is the reason we 



4 Strictly, this is an improper Riemann integral. 

5 We see Urn sup T _ 0+ dN/dr = 00 when dN/dr exists for 



< r. 



B. Divergence problem 

Let us consider the quantity tuf a which appears in the 
definition of Af [ flAip in Appendix]. We have 



-K iq E 3a 



e 2 N 



{iii - DiNj - DjNi)E 3 a . 



In order to examine this finiteness, we prepare 
1 (dt{e 2 ) 



Qa pi 



N 
1 

N 



d t {e 2 ) d 3 (e 2 ) 



^-^A J ~2d l {N l ) 



(3.2) 



Now we assume the conditions (a)-(d). Then N is finite, 

so the term 1/N is bounded below. When a parameter s 

is taken such that d s — {1, —N 1 } and s — corresponds 
the degenerate poi nt, the first and second terms in the 
right-hand side of (3^) are rewritten as, 



which diverges at s = 0. The di(N l ) term in (3.2) is 
finite since we assume condition (b). To sum up, we see 
( |3.2D diverges at the de gener at e po int. If we assume the 
triad reality conditions ( Al6 )-( A17 ), this fact tells us the 
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passing is impossible. Here the triad reality conditions 
are useful and stronger than the usual metric reality con- 
ditions as in Appendix A. 2. Since the triad is real, w° a 
is a real part of Af, so it is finite. And, since E z a is finite, 
ui® a E'a must be finite. The above conclusion contradicts 
this. Thus, unfortunately, we can not pass the degenerate 
point when the triad reality is assumed. 

Even if we do not impose the triad reality conditions, 
i.e. we impose the metric reality conditions or do not 
impose any reality conditions, this problem still exists 
as follows. In order to enforce the fmiteness of Af, its 
second term, ie a ( )c w l bc , must diverge such that it cancels 
out the divergence of ojf a . 

At the present stage, we do not know any exact so- 
lutions which satisfy this condition. Thus we conclude 
that the passing conditions (a)-(d), which aimed to pass 
a degenerate point directly, are difficult to satisfy simul- 
taneously due to this divergence of variable. 

IV. DEFORMED SLICE APPROACHES 

In the previous section, we showed that foliations 
through a degenerate point are difficult even within the 
Ashtekar framework. Therefore, in the second approach, 
we try to make a spacetime foliation avoid a degenerate 
point by breaking reality conditions. That is, we im- 
pose the foliations of 3-space which detour into complex 
spacetime only in the vicinity of such a degenerate point. 
This means that locally we are out of Einstein's frame- 
work. In this section, we examine what conditions are 
needed when we impose such a locally deforming time 
evolution. 

As we mentioned in the introduction, our aim is to keep 
continuous time evolution in the presence of a degenerate 
point. We extend the spacetime using complex numbers, 
and impose reality on the coordinates (t, x), so the folia- 
tion maintains 3+1 dimensions R 3 x R in C 4 . Note that 
our proposal of real coordinate is different from the Wick 
rotation, rather it is close to the recent proposal of com- 
plex lapse function by Hayward |Q. We use a freedom 
of gauge function to foliate into the complex spacetime. 

We assume a spacetime has a degenerate point and 
has also its future (a kind of coordinate singularity or 
focusing). In the following expression, we postulate that 
a single degenerate point exists at (t,x) — (t*,a;*), but 
extending to many degenerate points is straightforward. 
In order to recover the Einstein spacetime, the foliations 
have to satisfy the following conditions, such as the 'fo- 
liation recovering' condition which ensures that the foli- 
ation locally deforms into the imaginary region only in 
the vicinity of the degenerate point. 



A. Foliation recovering condition 

We seek the foliation which starts from the real sec- 
tion, evolves in the complex region only in the vicinity of 
a degenerate point, and ends in the real section. We call 
this the 'foliation recovering'. In the local deformation, 
the lapse function and the shift vector become complex 
valued to foliate in the complex region. Thus the 'foli- 
ation recovering condition' is expressed in terms of the 
gauge functions as: 

J J %N{t,x)dt = 0, (4.1) 
J + SN i (t,x)dt = 0, (4.2) 

where t— < t* and t + > t* are defined in appropriate far 
time from a degenerate point at t = t*. The reason for 
not imposing the other gauge function, the triad lapse 
Aq, is that the actual spacetime foliation is determined 
only by the lapse and the shift vector, and the triad lapse 
does not contribute in physical (external spacetime) foli- 
ations. 



B. Asymptotic reality condition 

Although we are considering deforming Einstein space- 
time by using a complex metric in the time evolution, 
we require that such deformations are as local as possi- 
ble. This is because our proposal is just a technique for 
treating a degenerate point. Therefore we impose real- 
ity conditions both at the spatial far limit and time far 
limit from a degenerate point, in order to ensure that 
the metric becomes complex only in the vicinity of the 
degenerate points. 

In the spatial far region from a degenerate point, the 
spacetime must be in real Einstein spacetime. In the 
time far limit, we also impose a condition to ensure that 
the metric becomes complex only in the time vicinity of 
the degenerate point. In other words, the initial value of 
E l a (t,x) is chosen to be real. Then we make time evolu- 
tion within Einstein spacetime before we meet a degen- 
erate point, and in the vicinity of a degenerate point we 
make time evolution into complex values, and return to 
the real values at the end. 

Those are expressed by gauge functions and metric 
variables. The former are: 

%N{t,x)^0 (4.3) 
SN%x)^>0 (4.4) 

for all four limits x — ► ± Aa;, t — * t* ± At where Aa; 
and At determine where we call the far region from a de- 
generate point, which depends on the particular problem. 
The latter is: 
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E l a E ja /detE{t,x) 



for all four limits x 



± Ax, t 



(4.5) 



U ± At. This 



expresses the primary reality condition ( All ), of which 
general discussions are in Appendix. 

As for the time-far limit of the metric reality condition 



(4.5), the expression can be written in an integrated form 



+ d 
~dt 



QEiE ja /detE(t, x)dt 



= 0. 



(4.6) 



The integrand is c alculated using the dyna mica l evolu- 
tion equation (AE), so the explicit form of (4.6) is not 
given at an arbitrary time except in some simple ana- 
lytic cases []. Normally, the condition ( |4.6| ) becomes a 
boundary problem of the evolution both in initial and 
final time-slice, so that any general application of this 
method requires numerical iterations for finding gauge 
functions to satisfy all conditions (4.1)-( |4~5| ). 

In general, we cannot deny t he p ossibility that two dif- 
ferent gauges satisfying (4.1 )-(4.5) exist and each yields 
a different real metric after deforming. However, there is 
no problem when we know the metric beyond the degen- 
erate point. In the next section, we show an example of 
such a foliation, where we fix our background metric to 
be flat spacetime with one degenerate point. 



V. EXAMPLES 

In this section, we discuss dynamical evolution through 
a degenerate point, using an analytic solution with the 
pull-back technique and a numerical demonstration. 



A. Exact solution with degenerate point 

One method to transform an exact solution is by the 
pull-back technique. When pull-back includes a singu- 
larity, it is not a diffeomorphism, and the coordinate 
transformation from the non-degenerate metric to the 
degenerate metric is singular. This technique was also 
used by Horowitz ||, when he discussed a possibility of 
topology change of the spacetime. Here, we derive a fiat 
Minkowski spacetime which has one degenerate point at 
the origin by using the pull-back method for later nu- 
merical demonstration of deforming slices. From the flat 
Minkowski metric, the following degenerate metric is con- 
structed: 



ds 2 = -[1 - (f(t)h(x)) 2 ]dt 2 

+ [-2f'(t)h(x)(l - f{t)h'(x))]dtdx 

+ (1 - f{t)h'(x)) 2 dx 2 + dy 2 + dz 2 , (5.1) 

where f(t), h{x) are an arbitrary functions which satisfy 

/(±oo) - /'(±oo) = 0, -1< f{t) < 1, 

and - 1 < f(t) < 1, (5.2) 

where the equality in the last equation is satisfied only 
when t = 0, and 



h(±oo) = h'{±oo) = 0, 
and - 1 < h'(x) < 1, 



■1 < h(x) < 1, 



(5.3) 



where the equality in the last equation is satisfied only 
when x = 0. 

We will see that this metric has the following proper- 
ties: (1°) flat, (2°) asymptotic Minkowski, (3°) degen- 
erate at t = x = 0, and (4°) Lorentzian (i-time, x, y, z- 
space) at elsewhere. 
Proof ) 

(1°) Let us put a = t, b = x — f(t)h(x). Then we see 
ds 2 = -da 2 + db 2 + dy 2 + dz 2 . 

(2°) Since /(*) -> and /'(<) -> when t -► ±oo, 
h(x) — > and h'(x) — ► when x — > ±oo, we see that the 
metric becomes flat for all four limits t — > ±oo, x — > ±oo. 
(3°) The determinant of this metric is detg^ = — (1 — 
f{t)h'(x)) 2 . Noting -1 < /(*) < 1 (equality is satisfied 
only when t = 0) and — 1 < h'(x) < 1 (equality is satis- 
fied only when x = 0), we obtain det g^ v < (equality is 
satisfied only when t = x = 0). 

(4°) We have g 00 = -1 + (f(t)h(x)) 2 . Since |/'(i)| < 1, 
\h(x)\ < 1, we obtain g 00 < 0. We easily see gu > 0. 
And we have det<? M „ < 0. Thus we see the metric is 
Lorentzian (t-time, x, y, z-space) except for t = x = 0. 

For exa mple, f(t) — e _t , h(x) = xe~ x satisfy condi- 
tions (5.2) and (5.3). We start from the metric given by 
the substitution of f(t) — e _t , h(x) = xe~ x into (5.1), 



ds 2 = -[1 - (2ta 
+4tae~* 2 ~ 



"* - x ) 2 }dt 2 

! [1 - (1 - 2x 2 ) e - 



^[1 - (1 - 2x 2 )e 



l dx 2 



~ x jdtdx 
dy 2 + dz 2 . 



(5.4) 



Hereafter, we call this the dM solution (a degenerate ex- 
pression of Minkowski spacetime). This metric is plane 
symmetric so that a degenerate point t = x = is not a 
point in spacetime. But this metric will be sufficient to 
show our proposal works in the next demonstration. 



Of course the reality conditions at an arbitrary time include 
this asymptotic condition. But reality conditions at an arbi- 
trary time are too restrictive to solve in the deformed slice 
approach. 



B. Numerical demonstration of two approaches 

Using the dM solution ( |5.4| ) as a background metric, 
we here show some numerical demonstration of our dis- 
cussions in the previous sections. The following com- 
putations were made by conventional finite difference 
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scheme using independent programs written in Fortran 
and Mathematica, both of which passed the appropriate 
convergence tests of grid size using the final imaginary 
part of density e. 



Since the metric (5.4) has a degenerate point at x — 
t = 0, it is clear that we can not evolve the solution 
beyond a degenerate point because the inverse 3-metric 
diverges at the degenerate point within the ADM for- 
mulation. As we denoted in §111, even if we take the 
Ashtekar formulation in time evolutions, some variables 
diverge at a degenerate point because of the 'divergence 
problem'. Thus we may say it is hard to 'pass' such a 
point. 

Fig. [l| shows such a feature. The figure is given by 

a time evolution of a solution (5T) using Ashtekar's dy- 
namical equation, and the real part of density e is plot- 
ted versus spatial and time comoving coordinates x and 
t. We set the gauge conditi ons (slicing conditions) the 
same as the exact solution (5.4), i.e., the lapse N = 1, 



the shift vector Ni = 2txe 1 
and the triad lapse Aq — 0. Initial data are constructed 



'[I - (1 - 2x 2 )t 



from the exact solution (5^4) taking triad e° as the diago- 
nal matrix (y/j xx , 1,1). We see that the time evolutions 
do not work properly after intersecting t = x = 0, be- 
cause of the 'divergence problem'. 




2.0 



Re (e) 



FIG. 1. An example of the intersecting approach for the 
dM metric (5.4) using Ashtekar's dynamical equations in evo- 
lutions. The density is plotted versus spatial and time co- 
moving coordinates x and t. A degenerate point exists at the 
origin (t = x = 0), and we see that the time evolutions do not 
work properly after intersecting the point. 

Next, we show that a 'deformed slicing' approach 
works well. In order to make a deformed foliation, we al- 
low the shift vector to be complex. The imaginary part s 
of N i ar e arbitrary but should satisfy conditions ( |4.2| ) 
and (4.4). As an example, we choose 



where a and b are arbitrary constants. We take the lapse 
N = 1 and the triad lapse Aq = 0, the simplest ones as 
before. 




Re (e) 
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Im(e) 



FIG. 2. An exam ple of the deformed slicing approach for 
the dM metric (5.4) using Ashtekar's dynamical equation in 



(5.5) 



evolution. The density [real part(a) and imaginary part(b)] is 
plotted versus coordinates x and t. A degenerate point exists 
at the origin (t = x = 0), and we see that the time evolution 
works properly for all the evolution and that the 3-metric 
satisfies the asymptotic reality condition again. 

In order to judge whether an evolution recovers asymp- 
totic flatness and satisfies asymptotic reality conditions, 
we checked two values at the final slice of our evolution; 

F(t final ) := max|5R(e(t = t finah x) - 1) | (5.6) 
R(tfinal) ■= max|3?(e(t = t finah x)) | (5.7) 

where these indicate the time asymptotic reality condi- 
tion (4.5). (F and R were used in convergence tests of 
the programs.) 

By changing a deformation of slice (two parameters a 
and 6), we searched for a solution which minimize the 
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flatness F and asymptotical reality R at the final time 
slice. Fig.|| shows a successful example of a local vio- 
lation of reality, in which we used a = 0.003982 and 



b = 2.5 in the imaginary part of shift vector (5.5), and 
set tinitiai = ~4,£ final = +4 in coordinate time. In the 
case of Fig.||, the flatness F and asymptotic reality R 
are F < 0(1O~ 2 ), R < O(10 -2 ) respectively when we 
take At = 0.01. We also should mention that this exam- 
ple satisfies asymptotic reality conditions at spatial far 
limit and foliation recovering conditions in the previous 
section. In the figure, density e is plotted versus coor- 
dinates x and t. We see that the time evolution does 
work properly in the sense that the real part recovers 
the analytic evolutions and the imaginary part vanished 
asymptotically. 

We found such a solution appears discretely in the 
two-parameter a-b plane. This eigensystem-like behav- 
ior is what we expected, since our requirements are over 
the freedom of this dynamical system as we discussed in 
§IV . This suggests that the spacetime dynamics in the 



complexified manifold will not necessary develop (or con- 
verge) into the real section of the manifold, without im- 
posing some reality recovering conditions. And note that 
a solution which satisfies asymptotic reality also satisfies 
asymptotic flatness. 

We also tried with other choices of shift vector, e.g. 
^N. t = atxe- b ^ t2 + x2 \ 3A 4 = 2atxe-^+^{\ - 2(1 - 
2x 2 )e~ b ( t +x )} and other forms where a and b are arbi- 
trary constants. We found that the general behavior of 
the discrete appearance of solutions is quite similar to the 
above example. The existence of solutions in the case of 
small b suggests to us that our deformed slicing approach 
is applicable not only for degenerate 'point' but also for 
degenerate 'region'. 



VI. DISCUSSION 

We have studied one of the advantages of Ashtekar's 
formulation: a tractability of degenerate points. We 
showed that a direct passing of degenerate point (which 
we call the intersecting slice approach) is hard to work 
in general if we impose some natural assumptions on dy- 
namical variables to express Lorentzian spacetime evo- 
lution. This means we can not make an evolution of 
hypersurface directly through a degenerate point even if 
we use Ashtekar's variables. 

Therefore we propose a deformed slice approach, which 
violates reality conditions locally near a degenerate point, 
and showed that this proposal works in an actual dy- 
namical problem. Therefore we conclude that a trick for 
dynamically passing a degenerate point exists in using 
complexified spacetime. This deformed slice approach is 
a natural foliation within Ashtekar's formulation since 
variables are originally defined as complex variables. 

Readers may think that if the essential trick of passing 
a degenerate point exists in complexifying variables, then 



we might find an example even within the complexified 
ADM formulation. We found also that this statement is 
true, and found similar discrete behavior in the deformed 
parameters. 

This result suggests to us that there is a new possibil- 
ity for dynamical problems in classical spacetime such as 
focusing or shell-crossing of coordinate points and also 
for singularities. Although this foliation requires param- 
eter tuning to satisfy reality recovering conditions, the 
fact that a real part of the solution always expresses the 
expected analytic solution indicates this dynamical evo- 
lution technique works for more general problems. 

There are some points which remain unclear when we 
break reality conditions locally. One big problem is a 
causality of this spacetime. Matschull |f|] commented 
on causality after a system passes a degenerate point. 
He classified the degenerate tetrad by rank, and consid- 
ers causal structure of each degenerate space-time in a 
real manifold. However we do not know whether such a 
usual causal structure can be extended to the complex 
manifold. 

If general relativity is extended to allow degenerate 
metrics, the topology of spacetime has a possibility to 
change even within a classical picture ||. If we allow 
breakings of reality conditions locally, then our classical 
path may connect the spaces with different topology dy- 
namically. We showed only the case that a spacetime 
will recover its original metric, but in general initial and 
final metrics are allowed to be different. Further analysis 
will broaden our research in dynamics which includes a 
degenerate point. 

A discrete appearance of solutions suggests the exis- 
tence of something like 'topological charge' of this sys- 
tem. We are now seeking such a charge, and plan to con- 
nect our discussions with a classically topology changing 
scenario and/or signature changing scenario in the early 
history of the Universe. 

We believe that this work is the first example to display 
some classical dynamical behavior of Ashtekar's variables 
for numerical evolution of spacetime. We expect our ap- 
proach will become the first step towards understand- 
ing dynamics of the signature changing process, topol- 
ogy changing process and causal structure in a complex 
manifold. 
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APPENDIX A: ASHTEKAR'S FORMULATION 

Here, we summarize the basic variables and equations 
in Ashtekar's formulation of general relativity, used in 
this paper. 



detE = ~e abc e ijk EiEiE*, where e ijk := e abc efe b e c k and 

, _— 1, |8| 

^ijk ' — e ^ijk |_|- 

Varying the action with respect to the non-dynamical 
variables N, N z and Aq yields the constraint equations, 



1. The Ashtekar variables 

The key feature of Ashtekar's formulation of general 
relativity Q is the introduction of a self-dual connection 
as one of the basic dynamical variables. Let us write the 
metric using the tetrad, e^, and define its inverse, 

E i, by g^y = fatvu and E i ztg^ilu- We define a 
50(3, C) self-dual connection 



± Aa ._ Oa _ ' a be 



(Al) 



where uij/ is a spin connection 1-form (Ricci connection) , 
u>* := E Iv V ^e^. Ashtekar's plan is to use only a self- 
dual part of the connection + A® and to use its spatial 
part + Af as a dynamical variable. Hereafter, we simply 
denote + Af. as Af. . 

Note that the extrinsic curvature, Kij = — (<V + 
n.;n')V;nj in the ADM formalism, where V is a covari- 
ant derivative on E, satisfies the relation —K^E^ 1 = oj® a , 
when the gauge condition E® — is fixed. So Af is also 
expressed by 



Af = -K^ a -e a bc u; bc 



(A2) 



The lapse function, N, and shift vector, N l , are ex- 
pressed as Eg = (l/Nj-W/N). This allows us to think 
of Eq as a normal vector field to E spanned by the con- 
dition t — x° =const., which plays the same role as that 
of ADM. Ashtekar treated the set (Af, El) as basic dy- 
namical variables, where E l a is an inverse of the densitized 
triad defined by 



El := eEl 



(A3) 



where e := det ef is a density. This pair forms the canon- 
ical set. 

In the case of pure gravitational spacetime, the Hilbert 
action takes the form 



S = I d 4 x[AfEl + -NElElF^ ab c - AN det E 



-WFZEl + A^E'l], 



(A4) 



where N := e 1 N, A is the cosmological constant [], 



V t El := d t El 



ie ab c A b E l c , and detE is defined to be 



C H = -e ab c ElElF^- Add E^0 } 

Cmi = -FfjEl w 0, 
Cca = ViEl 0, 



(A5) 

(A6) 
(A7) 



where F* v := {dA a ) llv --e a bc {A f\A c )^ v is the curvature 
2-form. 

The equations of motion for the dynamical variables 
(Af and El) are 

M = -ie ab c NE J b F?j + A" Fft + V t A^ + eANef , ( A8) 

Ei = -iD^NEiEi) + 2V j {N^Ej) + iA b eJEl, (A9) 

where 7V/};< : 0,1,!' - ieJA^Tf, for T % J + T? = 0. _ 
In order to construct metric variables from the vari- 
ables (Af,El,N,N l ), we first prepare tetrad Ej as 

Eft = (1/eN, —N i /eN) and E% = (0, E l a /e). Using them, 
we obtain metric such that 



f v := E»E v jr) U ' . 



2. Reality conditions 



(A10) 



Notice that in general the metric ( A10 ) is not real. 
In order to recover the real metric, we must impose the 
reality conditions. 

To ensure the metric is real- valued, we need to impose 
two conditions; the primary is that the doubly densitized 
contravariant metric 7 y := e 2 7 y is real, 



S*(E*#°) = 0, 



(All) 



and the secondary condition is that the time derivative 
of 7 U is real, 



%{d t (El& a )} = 0. 



(A12) 



Using th e eq uations of motion for El (A9), the gaug e 
constraint (A7) and the primary reality cond ition (All), 
we can replace the secondary condition (A12) with a dif- 
ferent constraint 



W ij := lfr{e abc E*E { jv k E$) « 0, 



(A13) 



7 We changed a factor in front of cosmological constant from 



1, e 
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which fixes six components of A® and E l a . Moreover, in 
order to recover the original lapse function N := Ne, 

we demand Q(N/e) = 0, i.e. the density e be real and 
positive. This requires that e 2 be positive, i.e. 



detE > 0. 



The secondary condition of (A.14) 



$S[dt(detJ5)] = 0, 



(A14) 



(A15) 



is automatically satisfied (see pi ) . Therefore, in order to 



ensure that e is real, we only require (A14). 

Rather stronger reality conditions are sometimes useful 
in Ashtekar's formalism for recovering the real 3-metric 
and extrinsic curvature. These conditions are 
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[13] L. Bombelli and R. J. Torrence, Class. Quantum Grav. 

7, 1747 (1990) 
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[15] H-J. Matschull, Class. Quantum Grav. 13, 791 (1996). 



and Q(Ei) = 0, 



(A16) 
(A17) 



and we call them the "primary triad reality condition" 
and the "secondary triad reality condition" , respectively. 
Using the equations of motion of E % a , the ga uge c onstraint 
(A7), the metric re ality conditions (All), (A12) and the 
primary condition ( |A16| ), we see that ( A17 ) is equivalent 

to a 



d i {N)E la + ~e 



'elNE^djEl 



(A18) 



From this expression we see that the second triad re- 
ality condition restricts the three components of "triad 
lapse" vector Aq. Therefore flAl8| ) is not a restriction 
on the dynamical variables (Af and E^) but on the slic- 
ing, which we should impose on each hypersurface. Thus 
the second triad reality condition does not restrict the 
dynamical variables any further than the second metric 
condition does. 



[1] A. Ashtekar, Phys. Rev. Lett. 57, 2244 (1986); Phys. 
Rev. D36, 1587 (1987); Lectures on Non-Perturbative 
Canonical Gravity (Singapore, World Scientific, 1991). 

[2] H. Kodama, Phys. Rev. D42, 2548 (1990). 

[3] G. T. Horowitz, Class. Quantum Grav. 8, 587 (1991). 

[4] E. Witten, Nucl. Phys. B311, 46 (1988). 

[5] G. Yoneda and H. Shinkai, Class. Quantum Grav. 13, 
783 (1996). 

[6] D.C Salisbury, L.C. Shepley, A. Adams, D. Mann, 
L.Turvan and B. Turner, Class. Quantum Grav. 11, 2789 
(1994). 

[7] I. Bengtsson, Class. Quantum Grav. 8, 1847 (1991). 
[8] M. Varadarajan, Class. Quantum Grav. 8, L235 (1991). 
[9] I. Bengtsson, Gen. Rel. Grav. 25, 101 (1993). 



9 



